Gate-controlled guiding of electrons in graphene 
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Ballistic semiconductor structures have allowed 
the realization of optics-like phenomena in elec- 
tronics, including magnetic focusing [lj and lens- 
ing [2j. An extension that appears unique to 
graphene is to use both n and p carrier types to 
create electronic analogs of optical devices hav- 
ing both positive and negative indices of refrac- 
tion [3j. Here, we use gate-controlled density 
with both p and n carrier types to demonstrate 
the analog of the fiber-optic guiding in graphene 
[4]-[8]. Two basic effects are investigated: (1) 
bipolar p-n junction guiding, based on the prin- 
ciple of angle-selective transmission though the 
graphene p-n interface, and (2) unipolar fiber- 
optic guiding, using total internal reflection con- 
trolled by carrier density. Modulation of guid- 
ing efficiency through gating is demonstrated and 
compared to numerical simulations, which indi- 
cates that interface roughness limits guiding per- 
formance, with few-nanometer effective rough- 
ness extracted. The development of p-n and 
fiber-optic guiding in graphene may lead to elec- 
trically reconfigurable wiring in high-mobility de- 
vices. 

Graphene is a single layer hexagonal lattice of carbon 
atoms with a gapless, linear dispersion that leads to novel 
electronic properties [9j [10]. Carrier type and density 
can be controlled via gates, creating high- mobility, bipo- 
lar graphene electronics [11] . Electronic transport across 
an interface of holes (p) and electrons (n) - a p-n junc- 
tion - has been studied experimentally using a combi- 
nation of top/bottom electrostatic gates to create p and 
n regions [T2UT4] . and is now well understood theoreti- 
cally [IMS]. 

An intriguing possibility is to use both n and p carrier 
types in bipolar graphene structures to create electronic 
analogs of optical devices having both positive and neg- 
ative indices of refraction. For example, a symmetrically 
biased p-n junction has been shown theoretically to cre- 
ate a negative refractive index medium analogous to a 
Veselago lens [3 . Such a device is not possible using 
conventional two-dimensional electron gas (2DEG) sys- 
tems, and demonstrates a unique feature of a Dirac-like 
band structure in conjunction with the ability to electro- 
statically tune between n and p carrier types. There is 
considerable theoretical interest in using the Dirac-like 
properties of graphene to create novel optical devices in 



graphene [4]-[8j I20H24] , but experiments on these systems 
have not yet been reported. 

In this Letter, we demonstrate experimentally and nu- 
merically the graphene analog of a well-known optical 
device, the fiber optic. Three regimes of current guiding 
are identified in the experiment and simulations: (1) p-n 
junction guiding, based on the principle of angle-selective 
transmission, (2) the graphene fiber-optic analog, using 
total internal reflection and (3) an mixture of the two 
effects. A metric of the guiding efficiency is predicted 
in all three regimes. By varying the external parameters 
of gate voltage and magnetic field, guiding efficiencies in 
each regime are extracted from experiment, where trends 
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FIG. 1: Schematics of device and guiding in unipolar and 
bipolar regimes, a Schematic of top-gated electron guiding de- 
vice with four contacts (i, c, gl and g2) used to measure resis- 
tances Ru, and r{ c . Voltages on top gate, Vtg> an d back gate, 
Vbg i ( n °t shown) independently control carrier densities (includ- 
ing sign), which serve as effective indices of refraction. Graphene 
lattice orientation is schematic and is not controlled, b, Optical 
guiding (OPG) is based on reflection above a critical angle when 
the density in the channel (under the top gate) is higher than out- 
side the channel (controlled by back gate), similar to the operation 
of a fiber optic, c, Alternatively, p-n guiding (PNG) is based on 
exponentially suppressed transmission through a p-n interface at 
oblique angles of incidence. 
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in the parameters are observed and compared to numer- 
ical simulations. 

Photons and electrons exhibit analogous wave phenom- 
ena, reflected in the similarity of the Helmholtz equation 
describing electromagnetic wave propagation and the 
Schrodinger equations describing propagation for elec- 
trons [25j|26]. In graphene, the Fermi energy (e) plays 
the role of index of refraction in an optical medium [Si- 
ll], with the important feature that e can be modified 
via electrostatic gates. We note that the dependence 
of the effective index of refraction on density and gate 
voltage in graphene differs from the dependence in con- 
ventional 2DEGs. In graphene, the wave number, and 
hence the effective index, is proportional to e. In conven- 
tional 2DEGs, the wave number and effective index scale 
as e 1 / 2 [2]. More significantly, for graphene, e can either 
be positive (for electrons) or negative (for holes). 

For the device shown schematically in Fig. la, the ef- 
fective index of refraction under the top gate (Region 
1), ei, is controlled by the combined voltages on the top 
gate, Vtg? an d backgate, Vbg? while the effective index 
of refraction outside the top-gated region (Region 2), 62, 
is controlled only by Vbg- When |ei| > | | , the de- 
vice operates as an electronic fiber optic with critical an- 
gle 6 ) c =sin _1 (|e2/ei|), This effect, termed optical guiding 
(OPG), is shown schematically in Fig. lb. Modes propa- 
gating with > C will be totally internally reflected and 
therefore travel down the channel without leaking out of 
the fiber. 

Transmission across a graphene p-n interface decreases 
exponentially with angle from the interface normal [T5] . 
Therefore at grazing incidence, nearly all carriers im- 
pinging on the p-n interface are reflected, which leads 
to guiding of all but a small number of carriers. This 
guiding mechanism is termed p-n guiding (PNG) and is 
shown schematically in Fig. lc. Depending on values of 
Vtg an d Vbg, the mechanism responsible for guiding will 
be either OPG, p-n -junction guiding (PNG) or a combi- 
nation of the two (OPG/PNG). The carrier-density loca- 
tion of these three regions are shown as a function of the 
density underneath the top gate (m) and outside the top 
gate (712) in Fig. 2a. For unipolar devices (no p-n junc- 
tions, shaded blue), only OPG is present. For bipolar 
devices (shaded red), both PNG and OPG/PNG can oc- 
cur. OPG occurs if 1 61 1 > \c2\; PNG occurs if e\ * 62 < 0; 
OPG/PNG occurs if both conditions are satisfied. 

Quantum transport simulations are used to extract 
guiding efficiency as a function of gate voltages (see 
Supp. Info, for numerical methods). The simulated de- 
vice (Fig. la) has four contacts, i (injector), c (collector), 
gl, and g2 (electrical ground 1, 2) (Fig. la). Guiding ef- 
ficiency is defined as the fraction of current collected at 
c due to injection from z, Q = Ti c /Tu, where Tmn is 
the transmission probability from contact m to n, and 
%i = Tic + %gi + 7lg2- Note that Q is finite even for 
equal indices of refraction in regions 1 and 2, since % c is 



\e 1 *e 2 <o 

N \ PNG 

opg\ 

PNG N \ 


/' OPG 
/I81IH82I 


I81IH82I / 
OPG / 


OPG 

\ PNG 
PNG X x 
6 1 *8 2 <6\ 




1.0 



FIG. 2: Simulated optical and p-n guiding in gated 
graphene. a Guiding regimes, OPG (blue), PNG and OPG/PNG 
(pink), as a function of density under the top gate (ni) and outside 
the top gate (112). The midpoint of the diagram is zero density in 
regions 1 and 2. b-d, Simulations of the current density J in the 
three regimes for fixed effective index of refractions e\ and 62- The 
x and y-axis represent the lOOnmxlOOnm size scale of the device. 
The guiding efficiency is 7 =0.53, 0.03, and 0.24 (Q = 0.58, 0.18 
and 0.29) in the OPG, PNG and OPG/PNG regimes, respectively. 
Green lines in (b-d) are lines of constant current density. 



nonzero when e\ = €2 i.e. a uniformly-biased graphene 
sheet. We therefore define 7 as the difference between ft 
and its value for equal indices, 



7(^1^2) = fi(ei, e 2 ) - ft(ei, ei). 



(1) 



The parameter 7(61,62) then serves as an effective mea- 
sure of the guiding efficiency due to unequal indices of 
refraction for the graphene channel, independent of the 
source of guiding. The condition fi(ei, ei) = Q(— ei, — ei) 
follows from particle-hole symmetry. We take Fermi en- 
ergy to be equivalent to the effective index of refraction, 
and so equate e\ and e 2 with Fermi energies in regions 1 
and 2. 

Simulations for a lOOnmxlOOnm device (match- 
ing the experimental geometry) in the OPG, PNG 
and OPG/PNG regimes (Fig. 2b-d), yield guid- 
ing efficiencies ft(-0.3eV, -0.2 eV) = 0.58 for 
OPG, ^(0.2eV,-0.2eV) = 0.18 for PNG, and 
^(0.3eV,-0.2eV) = 0.29 for OPG/PNG. Us- 
ing ft(0.3eV,0.3eV)=0.05 (Fig. 2b and d) and 
ft(0.2eV,0.2eV)=0.15 (Fig. 2c) gives 7 = 0.53 for OPG, 
7 = 0.03 for PNG, and 7 = 0.24 for OPG/PNG. To 
obtain these values, simulations assumed a root-mean- 
square interface roughness of 4 nm in the p-n case 
and no roughness in the p-p interface. Adding rough- 
ness is neccessary to obtain the qualitative trend ob- 
served experimentally, ft OPG >ft OPG / PNG >ft PNG ; with- 
out interface roughness for the p-n case, numerics gave 

fi OPG <fi OPG/PNG, Xhe kr 

ger roughness in the p-n 
regime presumably reflecting the poor screening of dis- 
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order at the zero-density p-n junction [27] and consis- 
tent with the theoretical observation that large-angle- 
scattering modes are deteriorated by disorder [28 . For 
p-p interface, a larger Q than experiments is obtained 
when assuming an ideal p-p interface. It is therefore nec- 
essary to also add some roughness to the p-p interface to 
obtain quantitative agreement with experiment, as seen 
in Fig. 3. 

We next discuss the experimental realization of elec- 
tron guiding. Devices were made from mechanical ex- 
foliation of highly-oriented pyrolytic graphite. Metallic 
contacts (i, c, gl and g2) were patterned with electron- 
beam lithography (see Fig la), and the size of device was 
reduced to 100 nmx 100 nm using an O2/AT2 plasma etch, 
giving a device of dimensions comparable to its mean-free 
path, determined by transport. A functionalized, top- 
gate oxide was grown [13] and a top-gate electrode pat- 
terned using electron-beam lithography. Differential re- 
sistance R = dl/dV was measured using a standard lock- 
in technique at a temperature T=30K. Relatively high 
temperature and large densities were used to suppress 
Coloumb blockade fluctuations of the resistance present 
in the small device at low temperatures and densities. 
The field-effect, Drude mean free path ^ m f p = h/e 2 -a /hp, 
where kp — \/7r\n s \, for micron-sized graphene sheets is 
routinely found to be > 100 nm away from the charge- 
neutrality point. Here, however, it is difficult to estimate 
£ m f p in these submicron devices, as the charging energy 
will create deviations from the simple Drude model. Ex- 
tracting -£ m fp from a four- terminal measurement ^(Vtg) 
yields ~70nm. Given the limitations of the Drude ^ m f p , 
the device is at very least in the quasi-ballistic regime. 

Experimental guiding efficiency fl e xp — % C /Tu is de- 
termined from transport measurements . Transmission 
Tu is obtained in a "channel" geometry, with current Ii 
applied to z, and c, #1, and g2 grounded, from the resis- 
tance Ru = Vi/Ii, 



h/2e 2 

Rii 



(2) 



Transmission % c is measured in a "focusing" geome- 
try pp, with current Ii applied to z, voltage measured at 
c, and gl and g2 grounded. Using this geometry, two 
resistances are measured by monitoring the voltages at 
contacts z and c: = Vi/Ii and R{ c = V c /Ii. A calcu- 
lation of ft exp could be made by taking the ratio of these 
two resistances, however, the symmetry % c = T C i and 
Tu — Tec, ma y n °t hold in real devices. Accounting for 
deviation from above idealities, we average ft exp over the 
two configurations Q exp = %i/Tu and Vt exp = T ic /T C c to 
obtain a value for guiding in terms of resistances R^ and 
R{ c (see Supplementary information), 
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FIG. 3: Effects of gating and disorder on guiding efficiency. 

a, (inset): Experimental guiding efficiency £l e xp (color scale) as a 
function of top-gate voltage, Vtg? an d back gate voltage, Vbg, 
extracted from the resistances Ru, and , with three guiding 
regimes indicated (black dashed lines). Cut line (white dashed line) 
through OPG regime where experimental and numerical guiding 
efficiencies are compared in main figure, (main): Experimental 
guiding efficiency 7exp, along constant fermi energy (effective index 
of refraction) in region 1, e\ = 0.3 eV. The values of 7 ex p rises from 
at Vqg — — 50 V — the point where the density is a constant in the 
device — to ~0.20 at Vqg =0V. Numerical guiding efficiencies are 
plotted for interface roughnesses of (black circles), 1 nm (green 
squares) and 2nm (blue diamonds), b, 7 ex p extracted along Vqg — 
— 10 V (red crosses), shows the extracted guiding in the OPG, PNG 
and OPG/PNG regimes along with 7 from simulation with disorder 
of (black) and 2nm (blue) in the OPG regime and (black), 2 
(blue), and 4nm (purple) in the PNG and OPG/PNG regimes. 
Good agreement between experiment is observed for 2 nm in the 
OPG regime and 4nm for the PNG and OPG/PNG regimes. 



where the two terms in Eq. [3] are T C i/T cc and Ti C /Tu 
respectively. Ideally, the symmetry of the device would 
entails T C i = % c and Tu = T cc and Eq. [3] would reduce 
to Vt = Tc I Tu- However, this will not be the case in 
typical experimental condition due to disorder and con- 
tact misalignment. A quantitative estimation of the ex- 
perimental guiding efficiency is obtained by taking the 
average of Q exp = T c i/Tu and Q exp = Tic/T cc , producing 
Eq. [3j Since < ft, 
4-7? z/ >0. R{-R 

with experiment (see Supplementary information). 



exp < 1, one also has the inequalities 
< R{ c and R{ c < R^ , in agreement 
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Values for ft exp as a function of top and back gate volt- 
ages, extracted Ra, and R{ c , are shown in the inset 
of Fig. 3, along with boundaries of the guiding regimes 
based on the boundaries in Fig. 2a. ft ex p is maximal in 
the OPG and OPG/PNG regimes, and follows the pat- 

tern ^>n^/™S«™ G - 

Experimental and numerical results for guiding effi- 
ciency are compared in the OPG regime, where guiding 
is most efficient, along a cut at constant density (Fermi 
energy ei~ 0.3 eV) in Region 1, indicated by the white 
dashed line in the inset of Fig. 3. Experimental guiding 
efficiency 7 exp along this cut is obtained from ft exp by 
subtracting ft exp (0.3eV, 0.3eV)=0.26. Figure 3 shows 
7 exp as a function of Vbg in the OPG regime along with 
numerical results for e\ = 0.3 eV with numerical inter- 
face roughnesses of 0, 1, and 2nm. Numerical values of 
7 were computed from ft using ft(0.3 eV, 0.3eV)=0.05. 
For all values of interface roughness, 7 increases with in- 
creasing Vbg (or decreasing | | ) , and good quantitative 
agreement with experiment is found for a roughness of 
2 nm. The dependence of 7 on Vbg can be understood 
by analogy with optical fiber, where decreasing the re- 
fractive index in the fiber cladding | | leads to smaller 
critical angle. The inset of Fig. 3 shows that the trend 
of increasing O exp with back-gate voltage is more promi- 
nent in the OPG regime than in the PNG regime. This 
is consistent with expectation, as the effect of Vbg in the 
PNG regime is mostly to change the location of the p-n 
interface, where reflection occurs; transparency of the p- 
n interface itself is not strongly affected by Vbg, as it is 
in the OPG regime. 

7earp is also extracted at Vbg=-10V, showing guiding 
in the OPG, PNG and OPG/PNG regimes (red crosses 
in Fig. 3b). Values for the "equal-episilon" background 
subtraction were extracted experimentally for the OPG 
regime, while the PNG and OPG/PNG values were in- 
ferred from the extracted value using particle-hole sym- 
metry (i.e. Q(ei,6i) = Q(— ei,— ei)) as in Fig. 2. In the 
OPG regime, j exp falls roughly linearly as a fucntion of 
V TG from 0.18 to 0, while in the PNG and OPG/PNG 
regimes it rises from 0.01 to 0.16. For Vtg between ~7 V 
and 10 V, j exp in not define as it does not satisfy the 
conditions for guiding (see Fig. 2a). Numerical simula- 
tions for this value of Vbg are shown for the 3 regimes 
for 3 different amounts of disorder: (black), 2 (blue) 
and 4nm (purple). Here we see good agreement between 
^y exp and 7 in the OPG for 2nm of disorder, similar to 
Fig. 3a. However 2 nm of interface disorder is clearly off 
in the PNG and OPG/PNG regimes, and more disorder 
(4 nm) is needed to obtain agreement between experiment 
and simulations. This is consistent with the argument of 
poor screening at the p-n interface used in the PNG and 
OPG/PNG regimes in Fig. 2. 

Transparency across a graphene p-n junction decreases 
with applied perpendicular magnetic field (B), as dis- 
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FIG. 4: Magnetic field improves gate-defined guiding, a, 

(inset) R{ c as a function of Vtg an d B at Vbg=-10 V. (main) Plot 
of R{ c at B=0T and B=5T. The increase in B corresponds to an 
increase in R^ c and the ratio R^ c (bipolar) /R{ c (unipolar), indicating 
an magnetic-field enhancement of PNG contribution to r\ c . Inset: 
R? c (VtGjB) demonstrates enhancement of R{ c is evident for fields 
B >2T. b, Simulation of J in the OPG/PNG regime at B=5T. 7 
is enhanced from 0.24 (Fig. 2d) to 0.50 (Q is enhanced from 0.29 
to 0.55). 



cussed theoretically by Refs. [I5j [29] . and seen exper- 
imentally in Ref. [12]. The reduced transparency in 
a magnetic field increases guiding efficiency, as demon- 
strated experimentally in Fig. 4. When R{ c 2 ^> R^R^ — 
Ra) — the experimentally relevant case, see Fig. SI in 
the Supp. Info — ft e xp ~ R{ c - R{ c (Vtg) a t B=5T is com- 
pared to the zero-field value in Fig. 4a. At Vbg=-10V 
an increase of ~0.1 kft in the unipolar regime and ~1 kft 
in the bipolar regime as B is increase to 5 T is observed. 
The inset of Fig. 4a shows i?{ c (VrG,B), where an enhance- 
ment in resistance is apparent for B>2T. Since the ratio 
R{ c (bipolar) /R{ c (unipolar) increases with B, we ascribe 
the enhancement in R{ c (and, by inference, an increase in 
^exp) as a result of an increase in the PNG contribution 
to current guiding. Fig. 4b shows the simulated current 
density with B=5T, showing an improved guiding effi- 
ciency of 7 = 0.50 (fi = 0.55) from the B=0 value of 
7 = 0.24 (ft = 0.29) [see Fig. 2d for B=0 value]. 

In summary, we have investigated electron guiding in 
graphene by tuning the carrier type and density using lo- 
cal electrostatic fields to create the analogue of an optical 
fiber. Guiding efficiency was extracted in three regimes: 
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in the unipolar OPG regime, where the device is analo- 
gous to a fiber optic; in the bipolar PNG regime, where 
guiding occurs because of reflection at p-n interfaces; 
and the bipolar OPG/PNG regime, where both mech- 
anisms operate. We also demonstrated experimentally 
that guiding efficiency increases with an applied perpen- 
dicular magnetic field, consistent with numerical results. 
One would expect that pseudomagnetic fields created by 
strain engineering [30 of greater than 300 T [31] can also 
enhance guiding. Improvements to guiding efficiency will 
result from reducing interface disorder. Engineering of 
a collimated source with modes within the acceptance 
cone of the fiber would also improves guiding. With such 
improvements, this approach could perhaps lead to elec- 
trically reconfigurable wiring using graphene. 
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